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An experimental review on the influence of nonlocal electrodynamics in the 
vortex pinning properties of non-magnetic borocarbide superconductors is pre- 
sented. We show that the pinning force density F p exhibits a rich and complex 
anisotropic behavior that sharply contrast with the small mass anisotropy of 
these compounds. For magnetic fields H applied parallel to the crystallographic 
c-axis, the first order reorientation transition between two rhombic lattices man- 
ifests itself as a kink in F p (H). For H_Lc-axis, a much larger F p (H) and a 
slower relaxation rate is observed. In this field configuration, nonlocality in- 
duces a fourfold periodicity in F p when H is rotated within the square basal 
plane. Unlike the out-of-plane anisotropy, which persists for increasing impurity 
levels, the in-plane fourfold anisotropy can be strongly suppressed by reducing 
the electronic mean free path. This result unambiguously demonstrate that the 
in-plane anisotropy is a consequence of nonlocal effects. 

1 INTRODUCTION 

The large attention that the borocarbide family of intermetallics (RNi2B2C, 
where R = rare earth) has received during the last years is in part due to the wide 
variety of interesting phenomena that these materials exhibit. Special attention 
has been devoted to R = Tm, Er, Ho and Dy where superconductivity and 
antiferromagnetism coexist in a large portion of the phase diagram. The high 
superconducting transition temperatures T c and the broad variation of the ratio 
Tn/T c (where Tm is the Neel temperature, ranging from 1.5 K to 10 K) make 
this family particularly appropriate to explore that coexistence P |2 |3J 

* Laboratorium voor Vaste-Stoffysica en Magnetisme, K. U. Leuven, Celestijnenlaan 200 
D, B-3001 Leuven, Belgium. 

tOak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6061 and Department of 
Physics, University of Tennessee, Knoxville, Tennessee 37996-1200. 

' Superconductivity Technology Center, MS K763, Los Alamos National Laboratory, Los 
Alamos, NM 87545, USA 

§Ames Laboratory, U.S. DoE and Department of Physics and Astronomy, Ames, Iowa, 
USA 

^Ames Laboratory, U.S. DoE and Department of Physics and Astronomy, Ames, Iowa, 
USA 



Even in the simplest case of non-magnetic borocarbides (R = Y, Lu), many 
remarkable properties have been reported. For example, in the superconducting 
mixed state the vortex lattice symmetry exhibits a strong field dependence with 
structural phase transitions separating a variety of exotic rhombic and square 
flux line lattices fFLL) P 171 151 131 1TU1 ITT1 IT21 1T3] . In addition, the equilibrium 
magnetization M deviates from the local London prediction^!] and exhibits 
a fourfold anisotropy when the field is rotated within the basal plane of these 
tetragonal materials [T51 1TC| . 

The presence of these non-hexagonal lattices has been attributed to the ef- 
fects of nonlocal electrodynamics, which arise when the electronic mean free path 
£ is larger than the BCS zero temperature superconducting coherence length £o- 
Nonlocal electrodynamics in superconductors had been traditionally associated 
with very low values of the Ginzburg-Landau parameter, k ~ 1. Borocarbide 
superconductors have k values in the range of 10 to 20, which allows one to 
avoid cumbersome core treatments by using the very simple London approach 
and simultaneously introduce the nonlocal effects as a perturbative term. 

This approach was first developed by Kogan et al. ^7J to analyze the case of 
a superconductor with an isotropic gap and an isotropic cylindrical Fermi sur- 
face. In this scenario, the authors successfully explained the observed deviations 
of the reversible magnetization M(H) of Bi:2212 in the mixed state from the 
logarithmic dependence, M oc ln(H c2 /H), predicted by the London model[T8|. 
Later on, Song et al.|Tl| found similar deviations in YNI2B2C when H || c-axis, 
that could also be quantitatively accounted for by this model. More recently 
we have extended that study to all field orientations and showed that this gen- 
eralization of the London theory provides a satisfactory complete description of 
the anisotropic M(H) with a self-consistent set of parameters^]. 

In general, the additional nonlocal corrections introduced in the London 
equation gives rise to a coupling of the microscopic supercurrents in a vortex 
to the underlying crystal symmetry. As a consequence, nonlocality induces an 
anisotropic response even in cubic materials [20]) which according to the local 
London theory, should behave isotropically. 

How actually the superconducting currents are connected to the crystal sym- 
metry is still a matter of debate. First, Kogan et al.(7JE] developed a model 
considering a superconductor with an anisotropic Fermi surface and an order 
parameter with s-wave symmetry. Later on, Franz and co-workers j2J analyzed 
the opposite case of an isotropic cylindrical Fermi surface in a d-wave supercon- 
ductor. The applicability of each one of these models depends eventually on the 
material studied. 

In the particular case of non-magnetic borocarbides compounds, there is 
growing evidence that both an anisotropic Fermi surface 22, 23. and a non con- 
ventional pairing mechanism [23 El ESI EHI 03 ESI ESI EOl EB are present. In the 
last years, the lack of any corroborative evidence for a non conventional pairing 
mechanism in these materials meant that most of the experimental results were 
interpreted within the Kogan et al.'s model. However, the proven existence of 
a pronounced anisotropy in the gap symmetry challenges the original view and 
strongly encourages one to critically re-examine several previous interpretations. 
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As we pointed out above, none of the introduced models has to be necessarily 
invoked in order to explain the deviation in the reversible magnetization from 
M oc ln(H C 2/H), whereas both of them correctly account for the anisotropic 
in-plane magnetization. However, still some differences appear between these 
approaches, when analyzing the field and temperature evolution of the vortex 
lattice. 

The model of Kogan et al. predicts |7] that two structural transitions in the 
FLL should occur in borocarbides for H || c: a first order reorientation transi- 
tion between two rhombic lattices at a field Hi and a second order transition 
from rhombic-to-square at H 2 > H\. According to this model, both transition 
fields shift up as temperature is increased. Recent small angle neutron scatter- 
ing (SANS) studies [TT1 1321 13"3") confirmed the existence of these transitions in 
YM2B2C and LUM2B2C. A jump in the apical angle (3 of the rhombic lattice, 
discontinuous within the resolution, occurs at Hi ~ 1 to 1.5 kOe, and the lattice 
becomes square (/3 = 90°) at H 2 ~ 2 to 2.5 kOe|35|. 

Franz et al. 21 showed that similar FLL rearrangements should appear in the 
case of a d-wave superconductor. However, two distinctive predictions emerge 
from this model. First, an additional first order transition should take place at 
very low temperatures. Second, in sharp contrast to the Kogan et al.'s model, 
the Hi boundary should decrease as temperature is increased. In a recent work, 
Nakai et al.|33| showed that when both anisotropic gap and Fermi surface are 
considered simultaneously, the intrinsic competition between these two effects 
leads to a much richer H — T phase diagram. 

Although the role of nonlocality on the equilibrium properties of the FLL 
has been convincingly established, much less attention has been paid to the 
effects on the dynamic or nonequilibrium vortex response. The very low critical 
current density J c observed in non magnetic borocarbides for H || c, associated 
with large pinning correlation volumes|8j, indicates that the elastic properties 
of the FLL must play a key role in the pinning. Since the shear modulus C 66 
depends on the vortex lattice symmetry [35], and the apex angle (3 undergoes 
a discontinuous jump at Hi, it is expected that Cqq and hence the pinning 
properties change abruptly at this fieldjHj. In other words, vortex pinning, 
which involves distortions from equilibrium vortex configurations, should be 
affected by the symmetry changes in the vortex lattice. Indeed, as we will show 
throughout this work, nonlocal effects influence the vortex pinning in several 
ways. 

First, in both non-magnetic compounds YNI2B2C and LUM2B2C, and for 
H || c, the reorientation transition at Hi induces a kink in the pinning force 
density F p (H). We find that H\(T) slightly decreases as T increases, in contrast 
to H 2 {T). We also study the effect of Co-doping in Lu(Ni 1 _ 2 ,Co ;r )2B2C. We 
observe that H\ decreases as the nonlocal effects progressively fade out when 
x increases. These results are in agreement with the expected evolution of Hi 
within a non conventional superconducting gap scenario. 

Second, we observe a fourfold oscillation in F P (H) when H is rotated within 
the ab- plane (in-plane anisotropy). We show that this effect is strongly sup- 
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pressed as the amount of impurities is increased, i.e. as the nonlocal effects are 
reduced. This clearly demonstrate that the in-plane anisotropy is a consequence 
of nonlocal effects. 

Strikingly, we also observe a large anisotropy in F p between the c-axis and 
the basal plane (out-of-plane anisotropy). We found that F p for H_Lc is about 
one order of magnitude larger than for H || c and has a quite different held 
dependence. Accordingly, a much faster normalized relaxation rate S is detected 
for H || c. This effect persists for all the explored impurity concentrations, 
demonstrating that the out-of-plane anisotropy is not related with nonlocal 
effects. We also rule out the presence of surface barriers for H_Lc by performing 
minor hysterisis loops. 

This work is organized as follows. In the next section we briefly introduce 
the experimental details and the samples studied. In section 3 we discuss the 
reorientational structural transition and how it influences the dynamic response. 
In section 4 we focus on the in-plane anisotropy of the pinning force density and 
in section 5 we discuss the out-of-plane anisotropy. Finally, in section 6 we state 
our conclusions. 

2 EXPERIMENTAL ASPECTS 

The composition, dimensions, critical temperature and estimated electronic 
mean free path for each of the studied single crystals is summarized in Ta- 
ble I. Normal state magnetization measurements in the Y-0 sample reflect the 
presence of a very diluted distribution of localized magnetic moments. The 
paramagnetic signal follows a Curie law which corresponds to a rare-earth im- 
purity content of 0.1 at. % relative to yttrium, probably due to contaminants 
in the yttrium starting material 1& . Similar measurements in the Lu-x single 
crystals show a much weaker Curie tail at low temperatures, T < 100 K, which 
might be due to a 0.001% magnetic impurities of Gd in the Lu crystallographic 
site 36 . 



Table 1: sample name, composition, volume, thickness, superconducting critical 
temperature and electronic mean free path for the investigated samples. 



sample 


composition 




V(mm 3 ) 


<(mm) 


Tc(K) 


1(A) 


Y-0 


YNi 2 B 2 C 




2.8 


0.5 


15.1 


300 


Lu-0 


LuNi 2 B 2 C 




2.5 


0.3 


15.7 


270 


Lu-1.5 


Lu(Ni 1 _ :c Co ;c ) 2 B 2 C (x= 


=1.5%) 


1.2 


0.4 


14.9 


100 


Lu-3 


Lu(Ni 1 _ a; Co ;c ) 2 B 2 C (x= 


=3%) 


0.44 


0.2 


14.1 


70 



Magnetization measurements were conducted on commercial Quantum De- 
sign SQUID magnetometers. Part of the relaxation measurements were taken 
in a MPMS-7 (H < 7 T) and the angular studies were performed in a MPMS-5S 
(H < 5 T) equipped with a home-made rotating system|53| IJSj. In the super- 
conducting mixed state, isothermal magnetization M(H) loops were recorded. 
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For the angular studies, both components of M(-ff) were measured using the 
two sets of pick up coils (longitudinal and transverse). The critical current 
density J c was determined from the width of the hysteresis loops according to 
the Bean critical state model, as previously described[3Sj- Relaxation measure- 
ments M(t) were performed on both branches (increasing and decreasing) of the 
hysteresis loop over periods of 1 hour. 

3 REORIENTATION AL TRANSITION 

According to the existing models, two types of structural transitions in 
the vortex lattice should occur in borocarbides materials, namely a first order 
reorientation transition between two rhombic lattices at a field Hi and a second 
order transition from the rhombic to a square lattice at a higher field H 2 . These 
predictions were recently confirmed by small angle neutron scattering (SANS) 
experiments |lll 1321 I33| on a sister crystal to the Y-0 sample studied here. For 
H || c and at T = 4.5 K, a jump in the apical angle (3 of the rhombic vortex 
lattice occurs at Hi ~ 1 to 1.5 kOe, and the lattice becomes square (j3 = 90°) 
at H 2 ~ 2 to 2.5 kOe. 

As we pointed out above, it is expected that the discontinuous jump of the 
apex angle (3 at Hi induces an abrupt change in the pinning properties through 
the shear modulus C66 |H1 ESI • I n order to check this prediction, we measured 
the pinning force density F p = |J C x B| for the Y-0 sample as a function of the 
applied field for H || c-axis at several T (see Figure^a)). 

We observe that at low fields F p (H) strongly decreases with increasing H, 
up to a certain field H* ~ 1.2 kOe, above which the field dependence is much 
less pronounced. Several facts indicate that this jump in dF p /dH is a signature 
of the reorientation phase transition Hi in the FLL0D]. First, the position of 
the kink for H || c coincides with the value of Hi reported in the literature[53j. 
Second, as shown in FigureEJ where F P (H) for sample Y-0 is plotted at several 
field orientations, H* is rather insensitive to the field orientation, in agreement 
with the behavior of Hi determined from SANS experiments |11|. 

In Figure ^a), it can be also observed that the transition field H*(T) re- 
mains almost constant at low T and, unlike the rhombic-to-square second order 
transition H2(T), it slightly decreases with T at higher temperatures. This be- 
havior is in agreement with recent neutron scattering experiments 33 and with 
theoretical predictions for a d-wave superconducting gap 21 . In contrast, it fol- 
lows from the model of Kogan et al. that both transitions, Hi and H 2 , should 
increase as nonlocality vanishes. In Figure|3we show the temperature evolution 
of H* together with the upper critical field H c2 . 

Further evidence that the observed feature at the field H* arises from non- 
local effects come from the dependency of this transition field with the elec- 
tronic mean free path. Indeed, the strength of the nonlocal perturbations is 
parametrized by a new characteristic distance, the nonlocality radius p(T,£): 
smaller p(T, £) corresponds to weaker nonlocality effects. This means that p(T, £) 
decreases with increasing T or decreasing I. Thus, if the position of the kink is 
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Figure 1: Field dependence of the pinning force density for H || c-axis for (a) Y-0 
sample at several temperatures and (b) Lu-0, Lu-1.5 and Lu-3 samples at T = 5 K. 

related with nonlocal effects, an increase in T or a decrease in i should produce 
qualitatively similar effects |71H7|. 

It has been shown that the Co substitutes on the Ni site in borocarbides 
compounds as a nonmagnetic impurity 36| that suppresses T c [37], but does not 
modify significantly the critical current J c 0]- This Co doping procedure allows 
one to accurately control the electronic mean free path by tuning the amount x 
of impurities. 

To determine the i dependence of H*, we performed measurements on the 
Lu-0, Lu-1.5 and Lu-3 samples for H || c (estimates of £ are given in Table I). 
The results are shown in Figure ^b), where we plotted F P (H) at T = 5 K. 
We observe that the kink in F p (indicated by the arrows) is still visible in the 
field range where the transition Hi should appear, and that it shifts to lower 
fields with increasing x, in agreement with the T dependence. In other words, 
the analogy between x and T is satisfied for H* , thus confirming that the kink 
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Figure 2: Field dependence of the pinning force density for three different angles of 
H with respect to the c-axis at T = 3 and 5 K. 

arises from nonlocality. 

From Figure H^b) we also confirmed the fact that the introduction of Co 
does not affect the pinning properties. Another interesting fact is that the Lu-0 
sample has a larger F p than the Y-0 sample at the same T, even though it has a 
lower density of magnetic impurities. This indicates that the magnetic moments 
in the Y-0 are not relevant pinning centers for the flux lines when H || c. 

4 IN-PLANE ANISOTROPY 

Let us now analyze the pinning properties for H_Lc. Figure 0] shows F p (H) 
for H applied along the symmetry axes [001], [100] and [110], at T = 5 K for 
the Y-0 sample. We observe that for the in-plane orientations, F p (H) exhibits 
a very different behavior than for H || c. First, over most of the held range 
F p (H.J-c) 3> F P (H || c). Second, the pinning force density exhibits a clear and 
systematic in-plane anisotropy, F p [100] > -F p [110]. We now focus on the origin 
of this in-plane anisotropy and defer the discussion of the large out-of-plane 
anisotropy to the next section. 

The fourfold nature of the in-plane anisotropy becomes evident in the lower 
panel of Figure where we show the angular dependence of the pinning force 
density, F p (ip), at T = 7 K for several H. This figure shows that for all H > 2 
kOe, F p (ip) exhibits a rather complex behavior: in addition to the main peaks 
at [100] and [010], secondary maxima are visible at [110] and [110]. 

For comparison, in the upper panel we show the in-plane angular dependence 
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Figure 3: Temperature dependence of the upper critical field and the structural re- 
orientation transition field for the Y-0 sample when H || c. 

of the Meissner slope —4n{dM/dH) for H = 30 Oe at T = 7 K. As H = 30 Oe 
is well below the lower critical field H c i for all ip, this curve represents the total 
flux exclusion of the Meissner state. The oscillatory behavior with periodicity 
7r (two-fold symmetry) originates from purely geometrical effects. Indeed, a 
field applied at any orientation within the basal plane can be decomposed in 
H x = Hcos(tp + 45°) and H y = Hsin(tp + 45°). If we approximate the crystal 
shape by the ellipse, the Meissner response associated with each component is 
4irMi = —Hi/ (I — Vi), where i = x;y and i>i = t/Li are the demagnetizing 
factors, thus 4?rM = -H[cos 2 ((p + 45°)/(l - v x ) + sin 2 (ip + 45°)/(l - v y )]. 

While it is very unlikely that any type of crystallographic defects could 
account for the complex in-plane variations in F p , nonlocal effects can provide 
a natural explanation for them. Due to nonlocality, the geometry of the vortex 
lattice depends on the orientation within the plane. We can again argue (as 
in the H || c case) that 6*66 depends on the lattice geometry, and that such 
dependence must be reflected in F p . Interestingly, F p has local maxima at 
the high symmetry crystallographic orientations [100] and [110], where Cqq is 
expected to exhibit local minima. In agreement with a nonlocal scenario, we 
have also observed that the in-plane anisotropy progessively decreases as T 
approaches T„[4Tj|. 

Conclusive evidence that the in-plane fourfold anisotropy arises from nonlo- 
cal effects is provided by similar measurements conducted in the doped samples 
Lu-1.5 and Lu-3. As we pointed out above, by tuning the impurity content in the 
sample, we are able to control the electronic mean free path and so the nonlocal 
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Figure 4: Pinning force density F p vs H with H applied along the principal axes 
[100], [110] and [001] at T = 5 K. 

effects. Figure IUb) shows the irreversible magnetization Af, ; (proportional to 
the pinning force density) as a function of tp for several fields at T = 5 K. Even 
though in this sample I is about three times smaller than in the Y-0 crystal, the 
fourfold anisotropy is still observed. For comparison, in Fig. Efa) we show the 
Meissner slope —AiT{dM/dH) = (1 — v)^ 1 which reflects the sample geometry 
effects. The smallness of the amplitude of this twofold oscillation indicates that 
the sample shape is closer to a disk. 

Interestingly, further decreasing the mean free path down to I = 70A by 
increasing the impurity concentration level to x = 3%, converts the oscillatory 
response from a fourfold to a twofold anisotropy (see Figure EJ) . Surprisingly, 
the irreversible magnetization does not maximize at the same angular position 
as the Meissner curve, but there is a shift of nearly 45° between them. This 
is an unexpected misalignment that cannot be accounted for by a bianisotropic 
current distribution, but may be related to the triangular shape of the crystal 
(viewed along the c-axis). 

5 OUT-OF-PLANE ANISOTROPY 

In the previous section we showed that the pinning force density in the Y-0 
sample is much larger for H_Lc than for H || c. This out-of-plane F p anisotropy 
sharply contrasts with the very small (< 10%) mass anisotropy [T51 IPol ITU) , 
therefore, explanations based on the anisotropic scaling frequently used in high 
T c superconductors@2ESI can be ruled out. 

In the same way, this effect is too large to be ascribed to nonlocality, which 
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Figure 5: In-plane angular dependence of (a) Meissner slope —A-kcLM / dH — (1 — v)~ X 
at T — 7 K (upper panel) and (b) pinning force density F p (ip) at T = 7 K and several 
fields (lower panel), for the Y-0 single crystal. The labels indicate the applied field in 
kOe. 



should appear as a perturbatively small effect. This conclusion is confirmed 
by the presence of similar anisotropy values in the doped samples Lu-1.5 and 
Lu-3, where nonlocal effects are strongly reduced. Indeed, in the main panel 
of Figure El we show F p (6), where 9 is the angle between the applied field and 
the c-axis, for the Lu-1.5 and Lu-3 at T = 5 K and H = 3 kOe. As the 
field is progressively tilted off from the c-axis, we first observe an almost angle 
independent pinning force F p up to 9 ~ 60°. Beyond this angle, F p increases 
very fast up to nearly one order of magnitud larger than for 9 = 0°, when H 
appoaches the basal plane. As we have previously demonstrated, this out-of- 
plane anisotropy is field dependent and drops to almost isotropic behavior near 
the peak effect (sec Fig.0J. 

One possible reason for the sharp increase of F p (0) as 9 — > 90° could be the 
presence of significant surface barriers for H_Lc. One way to check whether the 
surface barriers play a relevant role in the measured pinning force is to perform 
minor hysteresis loops with H || ab at several T and H. In the larger inset 
of Figure [5] we show an examples for the Y-0 sample at T = 5 K. If surface 
barriers were the main source of hysteresis, no flux changes would occur in the 
bulk while H is changing from one branch of the main loop to the other one, 
hence the data of the minor loop connecting the lower and upper branches would 
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Figure 6: In-plane angular dependence of (a) Meissner slope —A-kcLM / dH — (1 — 

at T — 5 K (upper panel) and (b) irreversible magnetization Mi(tp) at T = 5 K and 

several fields (lower panel), for the Lu-1.5 single crystal. 

be Meissner-like straight lines 03] • I n contrast, in the case of bulk pinning, the 
lines connecting both branches are curved (parabolic in the simplest Bean model 
for an infinite slab) just as we observe in the smaller inset of Figure^ Moreover, 
in a recent work|4*T] we have demonstrated that F p calculated from these minor 
loops assuming only bulk pinning arc in good agreement with those obtained 
from the main loops. Thus, a significant contribution to magnetic hysteresis 
arising from surface barriers can be ruled out. Fur her more, it is worth noting 
that the peak effect (near U c %) which is a bulk phenomena, is observed in 
the whole field orientation, clearly indicating that bulk pinning dominates the 
magnetic response |45|. 

A fact that we have not considered up to now is that, due to the significantly 
large time relaxation rates of the persistent currents in these materials, the cur- 
rent density J determined through magnetization measurements in the typical 
time scale of the SQUID-magnetometers (~ 20 sec) is smaller than the "true" 
critical current density J c . This suggests that the observed large out-of-plane 
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Figure 7: In-plane angular dependence of (a) Meissner slope —A-kcLM / dH — (1 — v)~ X 
at T — 5 K (upper panel) and (b) irreversible magnetization Mi(tp) at T = 5 K and 
several fields (lower panel), for the Lu-3 single crystal. 

anisotropy may be related to an angular dependence of the time relaxation of 
J. 

To explore this possibility, we have measured the normalized time relaxation 
rate of the irreversible magnetization S = —dhiJ/d\nt for the Y-0 and Lu-1.5 
samples for H || c and H±c at several fields. These measurements are shown in 
Figure[SJ Surprisingly, in both samples the creep rate S for H || c is comparable 
or even larger than the values obtained in HTSC at the same temperature (5 
K). In the simplest Anderson-Kim scenario, disregarding the possible influence 
of both quantum creep and glassy relaxation (see below), the pinning energy 
Up is given by U p rs T/S. For instance, in the particular case of the Lu-1.5 
sample at H = 6 kOe and T = 5 K, we obtain U p ~ 85 K, which is low even 
for HTSC standards. With regards to these values, it is worth keeping in mind 
that the scale of the elementary pinning energy u p is [H^/Bnn e£ 3 , where H c is 
the thermodynamic critical field and e is the mass anisotropy. In the HTSC, 
the low pinning energy is associated to the small u p , which in turn is due to the 
very small coherence length. In contrast, in the borocarbide crystals the small 
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Figure 8: Main panel: pinning force density F p as a function of the angle 9 between 
H and the c-axis, for the Lu-1.5 and Lu-3 samples at T = 5 K and H = 3 kG. The 
outer inset shows an hysterisis loop for H±c for the Y-0 sample at T — 5 K. A zoom 
in of the minor hysterisis loops corresponding to the full loop showed in the outer inset 
is shown in the inner inset. The dotted line shows the behavior expected for surface 
barriers. 



value of Up is a consequence of the very low density of defects. 

A remarkable fact observed in Figgis the large difference (almost one order 
of magnitude) of the creep rate between the two field orientations studied for 
the Lu-1.5 sample. We additionally find that for both orientations, and also for 
sample Y-0 (in the H \\ c orientation) S oc H n with an exponent n = 1.28(2). 
Considering again the simplest scenario, U p and F p are related by U p « F p V c £ 
where V c is the flux line lattice correlation volume. By combining this expression 
with S « T/U p , we obtain S « T/F p V c ^, so the creep rate anisotropy should 
be approximately the inverse of the F p anisotropy, as indeed observed. This 
result confirms that the source of the observed anisotropics is the anisotropic 
pinning energy U p . So, to complete the solution of this issue an explanation for 
the anisotropy in U p must be found. 

In brief, the origin of the large out-of-plane pinning anisotropy is unclear. 
A simple explanation for it could be the presence of some still unidentified 
anisotropic pinning centers, such as planar defects. Clearly, the complexity of 
the angular dependence of F p deserves further investigation. 

Another puzzling result comes from the study of the temperature dependence 
of the creep rate S(T) for fixed fields. In the inset of Figure OH we show S(T), 
for the Lu-1.5 sample at H = 4 kOe and H || c. We observe that S decreases 
monotonically as T decreases and it seems to saturate at low temperatures (T < 
4 K) in a value S ~ 0.02. Both the finite value of S as temperature approaches 
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Figure 9: Main panel: normalized time relaxation rate S(H) for different samples 
and field orientations. Inset: S vs. T for the Lu-1.5 sample at H — 4 KOe and H || c. 

to zero and its temperature independent behavior are typical indications of a 
quantum creep process. However, an estimation of normalized creep rate (near 
j c ) is given by Sq = fij (?7£ 2 L C ), T7: where -q is the friction coefficient for a vortex 
line and L c is the length of the tunneling vortex segment. This relation can be 
written as Sq ~ Q u yjc/jo where Q u = e 2 p n /he£,, p n is normal-state resistivity 
and e the mass anisotropy. For the Lu-1.5 sample, p n w 4.2 x 10~ 6 f2 cm, e ~ 1.1 
and jc/jo ~ 10~ 5 , therefore Sq ~ 4 x 10~ 6 , four orders of magnitude smaller 
than the measured value. The preceding analysis assumes that the dynamics 
of the quantum relaxation is dissipative. However, in our borocarbide crystals 
the electronic mean free path I is considerably larger than £o- This raises the 
possibility that the quantum relaxation may be more accurately described by 
a Hall- type vortex dynamics, which is supposed to dominate in the super-clean 
limit. In that case, an analogous expression for Sq should apply, with r) replaced 
by a Hall coefficient a 3> r?|47|. Thus, if that were the case, the numerical 
discrepancy with the experimental result would be even larger. 

An alternative interpretation could be that the observed T-independent be- 
havior corresponds to the well known plateau range predicted by the glassy 
relaxation and widely observed in HTSC. If that were the case, further reduc- 
ing the temperature (to T < 1 K) should result in a decreasing creep rate 
that tends to zero as T approaches zero. In this scenario, in the plateau range 
~ juln(l + t/to), and hence S ~ 0.02 implies that the glassy exponent 
/U w 2 (assuming, as is usually done, that ln(l + t/to) ~ 25). However, the 
p value obtained from the Maley analysis^Hj is only around 0.3. Moreover, a 
/i w 2 should produce very clear deviations from logarithmic relaxation, which 
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we have failed to observe. A natural way to define this issue would be to per- 
form a set of long term relaxation measurements (~ 12 to 24 hours). Another 
problem with this interpretation is that it implies an extremely small pinning 
energy, U p < T/iln(l + t/t ) for T > 2 K, that is, U p < 100 K. 

In summary, although the observed S(T) is suggestive of quantum creep, a 
quantitative description remains elusive at this point. The alternative expla- 
nation in terms of a glassy plateau also has serious difficulties. This complex 
situation highligths the relevance of the measurements of the normalized relax- 
ation rate S. Indeed, S is a fundamental variable of the dynamics of the vortex 
system, whereas the persistent current density J is measured at a particular 
time, and thus it is a more derived variable that depends on both the "initial" 
critical current J c and the relaxation rate. Clearly, more detailed studies of the 
flux creep phenomena in borocarbide superconductors are desirable. 

As a final application of the creep results, let us use them to estimate the 
dimensions of the correlation volume. As mentioned above, at H = 6 kOe and 
T = 5 K, we have U p ~ 85 K. At this field value F p ~ 2.12 kG.kA/cm 2 , which re- 
sults in a correlation volume V c « 10 10 A 3 . Because of the diluted distribution of 
pinning centers of these compounds, the Larkin-Ovchinnikov(LO) theory should 
apply without restrictions. According to this model V c = R 2 L C , where R c and 
L c are the dimensions of the vortex bundle perpendicular and parallel to the field 
direction and L c / R c ~ ^JlC^jC^, with C44 the tilt modulus and Cqq the shear 
modulus of the FLL. Within this description C44 « H 2 /4ir and in the simplest 
case of a triangular FLL in an isotropic medium with k 3> 1 in the intermediate 
field regime C§§ k H c \H{\ — h 2 )/l&-K^^. From this we can easily compute 
the size of the LO correlation length in the field direction L c ~ 10i? c ~ 10 4 A, 
which in units of the intervortex distance ao gives L c /ao ~ 20 for H = 6 kOe. 
This value turns out to be one order of magnitude smaller than that obtained 
by SANS measurements 8 at T = 2.2 K in Lu-0 and Y-0 crystals. However, we 
have naively assumed that C^s ~ H C \H /Air which according to the recent work 
of Knigavko et al. |49j is a overestimation of the total shear modulus. Indeed, in 
that work the authors demonstrate that the squash component Cqq decreases 
as H increases, vanishes at the transition field Hi and then grows slowly. In 
addition, Eskildsen et al. [HJ showed that for fields above 10 kOe the size of the 
correlation volume is dominated by the properties of the shear modulus. More- 
over, using the expression F p cx C^Cqq corresponding to the LO scenario and 
the experimental F P (H), we computed C§q{H), and surprisingly abtain that Cqq 
decreases monotonically with H both below and above Hi , in contradiction with 
the standard behavior |40j. This reinforce the idea that Cqq values derived from 
LO theory are larger than the "real" shear modulus and therefore the obtained 
L c tends to be somewhat smaller than those observed by SANS experiments. 



6 CONCLUSIONS 

The high quality of the borocarbides single crystals have the double benefit 
of producing large mean free paths and a very dilute distribution of vortex 
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pinning centers. The combination of these factors favours the observation of the 
influence of nonlocality on vortex pinning. In this work, we have shown that 
nonlocality significantly affects the irreversible properties of these materials. 
With the magnetic field directed along the c-axis, the pinning force density 
Fp(H) changes abruptly at the vortex reorientation field Hi, which is controlled 
by nonlocal effects. In addition, its dependence on temperature and electron 
scattering (mean free path) further ties the observed structure in F p to nonlocal 
electrodynamics. With the magnetic field in the square ct6-basal plane, the four- 
fold periodicity in F p and M irr (which is reminiscent of the four-fold periodicity 
in the equilibrium magnetization) again shows a coupling of the vortex lattice 
to the crystalline lattice; this higher order coupling (greater than a second rank 
mass-anisotropy tensor) arises from nonlocal effects, which again are washed 
out by thermal or electron-scattering disorder. These experimental observations 
demonstrate that nonlocality has a surprisingly profound influence on vortices, 
the flux line lattice, and vortex pinning in these remarkable materials. 

However, the rather complex angular dependence of the pinning force density 
clearly suggests that other mechanisms besides nonlocality are involved as well. 
In particular, the origin of the large out-of-plane anisotropy remains unclear. We 
have conclusively ruled out explanations based on the mass anisotropy, pinning 
by magnetic impurities, non-local effects and surface barriers. On the other 
hand, we have observed that this effect is also present in the time relaxation 
rate S, thus confirming that the source of the observed out-of-plane anisotropy 
is the anisotropic pinning energy U p . On top of this, the creep rate S for 
H || c exhibits two additional remarkable properties. First, the observed S 
values arc comparable or even larger than the values obtained in HTSC at the 
same temperature. Second, it is shown that S saturate at low temperatures 
in a finite value, likely owing to a quantum creep process. These unexpected 
behaviors cannot be quantitatively accounted for by the usual creep models for 
HTSC. Clearly, further experimental and theoretical studies are necessary to 
gain insight into this rather unexplored and fascinating field of the influence of 
nonlocality on the vortex pinning. 
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